Hydrodynamic afterburner for the Color Glass Condensate and the parton energy loss 
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We take hydrodynamic initial conditions in relativistic heavy ion collisions from the Color Glass 
Condensate (CGC) picture through the kr factorization formula. Gluon distributions produced from 
the CGC are found to provide good initial conditions for the hydrodynamic simulations in Au + An 
collisions at Relativistic Heavy Ion Collider (RHIC) energies. We reproduce the centrality, rapidity, 
and energy dependences of multiplicity within this approach. We also investigate the energy loss of 
high pt partons in the dense thermalized medium created from colliding two CGC's. We find that 
our results on the centrality dependence of nuclear modification factors for pions and back-to-back 
correlation for charged hadrons at midrapidity are consistent with the RHIC data up to semicentral 
events. Whereas, our approach in which jets are calculated from perturbative QCD 2^2 processes 
predicts less suppression at forward rapidity region compared to the BRAHMS data in Au-|-Au 
collisions at RHIC. 

PACS numbers: 



I. INTRODUCTION 

Experiments of high energy pp, dA, and AA colli- 
sions are currently performed at the Relativistic Heavy 
Ion Collider (RHIC) in Brookhaven National Laboratory 
for comprehensive understanding of the highly hot/dense 
matter, the quark gluon plasma (QGP). At collider ener- 
gies, collisions of two relativistic nuclei involve many as- 
pects of physics according to the relevant energy or time 
scale. There already exist many theoretical approaches 
to understand numerous RHIC data. We consider in this 
work, particularly, the physics of the gluon saturation 
in a colliding nucleus, hydrodynamic evolution of pro- 
duced matter, and the energy loss of hard partons in the 
medium. Our goals are to combine them and to take a 
first step toward a unified understanding of the dynami- 
cal aspect of high energy heavy ion collisions. 

First data reported by the STAR Collaboration at 
RHIC 0j revealed that the observed large magnitude of 
elliptic flow for hadrons is naturally explained by hydro- 
dynamics. This suggests that large pressure possibly in 
the partonic phase is built at the very early stage (r 0.6 
fm/c) in Au+Au collisions at RHIC. This is one of the 
strongest indications of early thermalization of the QGP 
at RHIC ^|3||. These hydrodynamic predictions also give 
good agreements with mass dependences of the second 
harmonic coefficient of azimuthal distribution V2{pt) in 
the low Pt region . This is in contrast to the case in the 
lower energy collisions where the hydrodynamics always 
overpredicts the data [Ml^' Hydrodynamics also predicts 
that the scaled elliptic flow, which is defined as the sec- 
ond harmonics V2 divided by initial spatial eccentricity e, 
becomes almost constant around 0.2 Q. The experimen- 
tal data reaches the hydrodynamic limit for the first time 
in central and semicentral collisions at RHIC energies @ . 
Moreover, V2 as a function of pseudorapidity [9| can be 
reproduced by hydrodynamics only in | 77 |< 1 All 
these analyses indicate that a high dense matter created 



in heavy ion collisions at RHIC energies achieves (local) 
thermal equilibrium state in the low pT region, around 
midrapidity, and up to semicentral collisions. 

On the other hand, jet production rate becomes large 
enough to measure high px power-law tails at RHIC. 
These minijets go across the expanding matter rather 
than participate in the dense medium in the case of heavy 
ion collisions. During traveling through the expanding 
medium, energetic partons interact with soft matter and 
lose their energies dynamically (jet quenching). Thus, 
high Pt hadrons coming from fragmentation of minijets 
can be good probes of the bulk matter It was ob- 
served that nuclear modification factors for hadrons in 
the high pT > 5 GeV/c region are considerably smaller 
than unity in central collisions 0, 0, Q, 0| and that 
it increases with decreasing centrality and eventually 
reaches around unity in peripheral collisions. Disap- 
pearance of the away-side peak in azimuthal correlation 
functions for high pT hadrons is also observed in cen- 
tral Au-I-Au collisions 0|. This can be described by 
perturbative QCD (pQCD) based parton models with 
parton energy loss in the dense medium 0, 0, 
Cronin enhancement of the hadron spectra and existence 
of the back-to-back correlation at mid-rapidity in dA col- 
lisions at RHIC [T5ll20l l2ll[2^ support that suppression 
of the yield and disappearance of the away-side peak 
in Au-|-Au collisions at mid-rapidity are due to the fi- 
nal state interaction. The radiative parton energy loss 
induced by the medium has been studied extensively 
m El il El m Ell and found that the energy loss is 
proportional to the gluon rapidity density of the medium. 
Theoretical analyses reveal that a large parton density is 
required to account for observed supra-essions of single as 
well as dihadron spectra at RHIC [l3.l2^ l30j. 

From the above discussion, the current RHIC data 
strongly suggest that the parton density created in heavy 
ion collisions at RHIC is dense enough to cause both 
large elliptic fiow of bulk matter and large suppression of 
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high-pT hadrons. One of the important problems in the 
physics of heavy ion colhsions is to reveal the origin of 
this dense matter. The bulk particle production in high 
energy hadronic collisions is dominated by the small x 
modes in the nuclear wave function, where x is a momen- 
tum fraction of the incident particles. Therefore, an ori- 
gin of the large density could be traced to the initial par- 
ton density at small x inside the ultra-rclativistic nuclei 
before collisions jsj ■ It is well known that gluon density 
increases rapidly with decreasing x by the BFKL cascade 
until gluons begin to overlap in phase space js^l where 
nonlinear interaction becomes important. These gluons 
eventually form the Color Glass Condensate (CGC) [33l |. 
This phenomenon is characterized by a "saturation scale" 
which corresponds to the gluon rapidity density per 
unit transverse area. When is large {Qs ^ Aqcd), 
the coupling becomes weak {as{Qs) ^ 1). In addition, 
the gluon occupation number in the wave function is 
huge, ~ l/a.s{Qs) ^ 1- Therefore this can be studied by 
an weak coupling classical method known as McLerran- 
Venugopalan (MV) model [sJl • It has been developed to 
compute the classical parton distributions of nuclei based 
on this model. Moreover, renormalization group methods 
which systematically incorporate quantum corrections to 
the classical effective theory were developed HE HE 113 ■ 
There is an evidence for gluon saturation in deep inelas- 
tic scattering at HERA. Models based on the saturation 
picture fit the HERA data well in the moderate and 
low X [sE 113 • This effective theory can be applied to the 
collisions of nuclei and its solution can be used for ini- 
tial conditions of further evolution of the system created 
in heavy ion collisions 0, 0, 0, 0| . Perturbative so- 
lutions of the gluon production for collisions of two nu- 
clei ^E QEl as well as pA collisions |46l |47| | have been 
computed. Non-perturbative classical solutions of the 
Yang-Mills equations were also obtained by the real-time 
numerical simulations on the lattice 

m, EES- The 

calculations based on the CGC picture have been com- 
pared to various RHIC data. Remarkably, the CGC re- 
sults on the global obscrvables, e.g., the centrality, ra- 
pidity and energy dependences of charged hadron mul- 
tiplicities agree with the RHIC data under the assump- 
tion of parton- hadron duality [HEl- The CGC picture 
is also applied for studies of particle correlations and az- 
imuthal anisotropics |53 | , and valence quark distributions 
at small- a; [sE ■ It has been shown that the classical wave 
function in the MV model contains Cronin enhancement 
and that quantum evolution with respect to x makes 
the spectrum suppressed [H^ . l7(i| | . Suppression of nuclear 
modification factor at forward rapidity region in dA col- 
lisions [IE in the CGC picture is also consistent with the 
recent RHIC data 

The above facts show that the CGC, hydrodynamics, 
and the energy loss of hard partons are the key ingre- 
dients to describe the RHIC physics. One expects that 
they are indeed closely related with each other. For ex- 
ample, the CGC could be a good initial condition for a 
thcrmalizcd state because it produces a large number of 



gluons. Thus a large number of gluons are responsible 
for large suppression of jet spectra. In order to under- 
stand systematically the whole of the stages of relativis- 
tic heavy ion collisions, it is indispensable to construct a 
unified and dynamical framework in which the above dif- 
ferent ideas are consistently incorporated. In this work, 
we perform hydrodynamic simulations with the initial 
condition taken from the CGC and simulate parton en- 
ergy loss in the fiuid elements. This approach will lead 
us to get deeper understanding of the dynamical aspect 
of heavy ion collisions. 

In addition, some of the problems which are inherent 
in a particular approach can be removed. For instance, 
one conventionally parametrizes initial conditions in the 
hydrodynamic simulations shortly after a collision of two 
nuclei so as to reproduce final observed spectra. Hydro- 
dynamic framework itself does not provide an initial con- 
dition. Hence one has to choose it from infinite possible 
candidates. Any initial conditions which reproduce the fi- 
nal hadron spectra are equivalent, so it is hard to discard 
even one of them within the hydrodynamic approach. 
This is the ambiguity problem of hydrodynamic model 
in heavy ion collisions. A good example can be found in 
Ref. |57l |: Two completely different initial conditions end 
up almost the same rapidity distribution. Therefore it is 
desired to take an initial condition which is obtained by 
a reliable theory. On the other hand, most of the cal- 
culations based on the CGC do not include final state 
interactions. According to estimation of initial parton 
productions from a classical lattice Yang-Mills simula- 
tion based on the CGC HE EE; the transverse energy per 
particle is roughly Et/N ^ 1.5 GeV, while the experi- 
mental data yields Et/N - 0.6 GeV HE- The elliptic 
flow from a classical Yang-Mills simulation is inconsis- 
tent with RHIC data [sE- Those discrepancies, which 
are mainly due to the lack of collective effects, will be 
removed by introducing further time evolutions through, 
e.g., hydrodynamics. 

For the calculation of parton energy loss, one needs 
time evolution of parton density. Bjorken expansion 
is assumed for the evolution of parton density in almost 
all work except Ref. (62j in which hydrodynamic simula- 
tions are used to estimate the energy loss of partons for 
the first time. We have previously established a model 
(the hydro-|-jet model) in which hydrodynamic evolutions 
are combined with explicitly traversing non-thermalized 
partons with large transverse momenta 0, [IE 0, EE ■ 
Here the production of high pt partons is described by a 
pQCD parton model. These high pT partons come from 
parton distribution function in relatively large x region. 
Systematic studies based on th e hy dro-f-jet model have 
been performed for pT spectra |65l |. back-to-back corre- 
lation functions pseudorapidity dependence of nu- 
clear modification factors , and elliptic flow 0, EE ■ 
For the initial condition of hydrodynamic simulations, we 
simply parametrized the initial shape of energy density 
in the previous work. Instead, we will employ the CGC 
picture and use gluons produced from melting CGC for 
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initial conditions for hydrodynamic simulations in this 
paper. The number and the energy density distribu- 
tions for produced gluons are calculated from the fc^ 
factorization formula [HJl- Although the produced glu- 
ons will reach a thermalized state through the dissipa- 
tive processes in the realistic situations, the description 
of non-equilibrium phenomena is beyond the scope of the 
present paper. Instead, we assume thermalization for 
gluons produced from the CGC. 

This paper is organized as follows. In Sec. ^ we 
discuss how the gluon distribution produced from the 
CGC is used for the initial condition in the hydrody- 
namics. In Sec. mil we briefly summarize the hydro-|-jet 
model. In Sec. IIVI we present hydrodynamic results on 
(pseudo)rapidity distributions for charged hadrons and 
the transverse momentum distributions for pions, kaons, 
and protons. We show that the CGC provides good ini- 
tial conditions in the hydrodynamic model. We also show 
that, with almost keeping the shape of rapidity distri- 
bution during the hydrodynamical evolution, transverse 
collective flow is generated by pressure gradient perpen- 
dicular to the collision axis. We study the centrality 
dependence of the nuclear modification factor for neu- 
tral pions and the back-to-back correlations for charged 
hadrons. The previous analysis on the high pT hadron 
spectra in the forward rapidity region is revisited. In the 
final section, we give a summary and a discussion about 
the further improvements in our approach. 



II. INITIAL CONDITION FOR 
HYDRODYNAMICS FROM GLUON 
SATURATION 



There exists a lot of extensive work for the description 
of gluon productions in nuclear collisions in the satura- 
tion regime where nonlinear effects become important. 
Perturbative solutions for the collision of two nuclei in 
the MV model were obtained in Refs. [3 [i^. How- 
ever, those solutions are limited for the relatively high 
momentum and cannot be used for the calculation of to- 
tal gluon multiplicities. In Ref. [g^, analytic expression 
for the gluon multiplicity is obtained. Non-perturbative 
classical boost invariant solutions are currently available 
on the lattice and it is highly desired to use 

the lattice results for obtaining reliable initial conditions. 
However, one of our purposes in this paper is to investi- 
gate the rapidity dependence of final hadrons. Therefore, 
we employ the kr factorization formula along the line of 
work done by Kharzeev, Levin, and Nardi (KLN) jsj for 
the computation of the gluon rapidity distribution which 
is to be the initial condition for sequential hydrodynamic 
evolution. 



A. Gluon Production in the kr factorization 
formula 

The number of produced gluo ns in the A:T-factorization 
formula is given by IstL Isslls^ 

d^y ^ i^y"-"^(^) 

X (j)Aixi,kT;X±_)(l)B{x2, (pt - kTf]Xj_), 

(1) 

where xi^2 = Pt ^^pi^y) / Vs and y and pr are a rapid- 
ity and a transverse momentum of a produced gluon. 
Running coupling ag is evaluated at the scale = 
max(A:|,, {pT — fcr)^)- The unintegrated gluon distribu- 
tion (j) is related to the gluon density of a nucleus by 

xGa{x,Q'^) = / d^kTd^x±(j)A{x,kT;x±). (2) 

There are several models/parametrizations in the liter- 
ature which provide a gluon distribution in the satura- 
tion region. In principle, the x dependence of uninte- 
grated distribution functions in Eq. should be solu- 
tions to the nonlinear quantum evolution equations such 
as Balitsky-Kovchegov equation iTCt l7ll | or, more gener- 
ally, JIMWLK equation I^E l37| . Leaving this analy- 
sis for the future work, we follow basically the KLN ap- 
proach which captures the essential features of the gluon 
saturation physics in a simple manner. 

In the MV model the gluon distribution below 
the saturation scale is logarithmically suppressed 
^(fc|,) ^ ln{Ql/k^) compared to the perturbative dis- 
tribution ~ Instead, motivated by the KLN ap- 
proach, we use a simplified assumption about the unin- 
tegrated gluon distribution function: 

=|^#' '^f^^' (3) 

where Cp = {N^ — l)/(2iVc). We introduce a small reg- 
ulator A = 0.2 GeV/c in order to have a smooth distri- 
bution in the forward rapidity region \y\ > 4.5 at RHIC. 
Other regions are not affected by introducing a small reg- 
ulator. As we will discuss below, the above distribution 
depends on the transverse coordinate x± through Q^. 
We are interested in the centrality and rapidity depen- 
dences of total yields which arc dominated by the low 
transverse momenta. The above simple form should be 
fine for this purpose, although the effect of anomalous di- 
mension resulting from the presence of the CGC modifies 
the Pt dependence of the produced gluons as shown in 
Ref. Ill El III El m 

An overall constant factor k is 
determined by fitting the multiplicity of charged hadron 
at midrapidity at y'sjvw = 200 GeV for the most cen- 
tral collision. Saturation momentum of a nucleus A in 
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A + B collisions is obtained by solving the following im- 
plicit equation with respect to Qg at fixed x and x± 



2n 



where 



p^,,d^^) = Ta{x^) {1 - [1 - a'j^Mx±)/B]^} . (5) 
The number of participants is given by 



part 



C«'a:^-L(Ppart(a:i)+Ppart(a'±))- (6) 



We take a simple pcrturbativc form for the gluon distri- 
bution for a nucleon 



xGix,Q^)= Kin 



^QCD 



(7) 



where A = Aqcd = 0.2 GeV. K is used to control sat- 
uration scale in Eq. (0} [t^. We choose K ~ 0.7 for 
A = 0.2 {K = 0.55 for A = 0.25 or if = 0.44 for A = 0.3) 
so that average saturation momentum in the transverse 
plane yields {Ql{x = 0.01)) ~ 2.0 GeV^/c^ in Au-t-Au 
collisions at impact parameter 6 = at RHIC. Similar 
to the KLN approach, x~^ dependence of the satura- 
tion scale is motivated by the Golec-Biernat-Wiisthoff 
model [S^- The factor (1 — x)" shows that gluon density 
becomes small at x — *■ 1. n usually depends on the scale 
Q^. Here we take n = 4 as in the KLN approach [Flj . 

We obtain the rapidity distribution for produced glu- 
ons at each transverse point x± by performing the inte- 
gral of Eq. numerically. The transverse energy distri- 
bution dEx/dy is also obtained by weighting the trans- 
verse momentum of gluons pr in the integration with 
respect to pt in Eq. ^ . We cut off the integral range of 
Px in Eq. since only the low px partons are assumed 
to reach the local thermal equilibrium. Instead, hard 
partons are included by assuming usual pQCD 2 — > 2 
processes in which DGLAP evolution is included as it is 
important in the high reactions. It is unclear which 
value we have to use for the cut off momentum pr, cut, 
because it contains a nonperturbative nature and also it 
would depend on the choice of unintegrated wave func- 
tion. We set the cut off momentum at pr.cut = 3 GeV/c 
which corresponds roughly to the maximum saturation 
scale at X = 0.01 at the origin x± = in central Au+Au 
collisions at RHIC. 

Figure ^(e^) shows rapidity distributions of produced 
gluons obtained from Eq. Assuming free stream- 

ing in the longitudinal direction until r = 0.6 fm/c, we 
also obtain the number density and the energy density 
as functions of a transverse coordinate as shown in Fig. Q] 
(b) . Here the direction of the a;-axis is the same as that of 
the impact parameter vector. Solid (dashed) line corre- 
sponds to the number (energy) distribution of produced 
gluons. 




FIG. 1: (a) Rapidity dependence of initial gluon transverse 
energy (dashed line) and number distribution (solid line) in 
Au + An collisions at ^^sWn = 200 GeV at 6 = 2.4 fm. (b) 
Gluon number and energy densities as functions of a trans- 
verse coordinate. Parameters are = 3.6, K = 0.7, and 
A = 0.2. As we will discuss in the next subsection, these 
parameters correspond to an initialization of hydrodynamic 
simulations which matches the number density of gluons pro- 
duced by collisions of two CGC's with the parton distribu- 
tion in hydrodynamic simulations at initial time. The value 
of is chosen so that we reproduce the multiplicity of final 
hadrons as discussed in Sec. IV. In this procedure, the trans- 
verse energy distributions represented by dashed lines are not 
used as hydrodynamic initial conditions. In the next sub- 
section, we will also discuss the other matching procedure in 
which the transverse energy distribution is used instead of the 
gluon number density. In that case, the gluon distributions 
are factor 1.6 smaller than the results in these figures. 



Regardless of the range of the parameter A = 0.2-0.3, 
we find kx factorization formula gives almost the same 
shape of rapidity distribution. The initial transverse en- 
ergy per particle is estimated to be Ex/Ng ~ 1.6 GeV at 
y = 0. This is within a range estimated in a numerical 
simulation of the classical Yang-Mills equation [sol . Is^ . 
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B. How to connect the CGC to hydrodynamics 

In order to obtain a correct initial condition for hydro- 
dynamics, one needs non-equilibrium description for the 
collisions of heavy nuclei. Especially, whether thermal- 
ization is achieved within a few fm/c is a longstanding 
problem in the physics of heavy ion collisions. The Boltz- 
mann equation with only the elastic parton-parton scat- 
tering predicts that the system never reaches a thermal 
state within a reasonable time scale in heavy ion colli- 
sions ^3 • Inelastic processes such as gg ggg are quit e 
responsible for the equilibration of the system |4ll. l4,3l| . 
Although the thermalization process can alter initial con- 
ditions, the centrality dependence of total hadron mul- 
tiplicities is approximately preserved in the CGC ini- 
tial conditions |42|. QCD plasma instabilities may also 
play an important role in the thermalization of matter 
from anisotropic initial conditions and these may help 
to speed up the equilibration of the plasma . A ques- 
tion of chemical equilibration has been addressed in Refs. 

[zilzi. 

Leaving these problems for the future work, we as- 
sume that the system of gluons initially produced from 
the CGC reaches a kinematically as well as chemically 
local equilibrate state at a short time scale. We further 
assume that during being thermalized, the shape of the 
rapidity distribution is not changed. Thus, we take ini- 
tial conditions from gluon distribution obtained in the 
previous subsection based on the CGC. 

We consider two simple ways to provide initial condi- 
tions from the CGC, i.e., matching of number density 
(IC-Ti) and matching of energy density (IC-e). We will 
compare charged particle distributions from these ini- 
tial conditions with each other in Sec. IIV Al Assuming 
Bjorken's ansatz y = rjs where 77s is the space-time 
rapidity 77s = (1/2) ln(a;~'"/a;~), we obtain the number 
density and the energy density for gluons at a space-time 
point (ro,a;^,?7s) = {to,x) from Eq. JQ) 



ng{To,x) : 
eg{To,x) 



Tod%d'^x± ' 
TQdi]sd'^xj_ ' 



(8) 
(9) 



We take two steps to specify initial conditions in hydro- 
dynamic simulations by using Eqs. |(SJ) or (O. Firstly, 
we assume that these partons are thermally equilibrated. 
Secondly, the gluon number density simply represents the 
Nf = 3 parton number density, i.e., the QGP. Although 
the initial time tq could depend on incident energies 
(or even on the space coordinate through the saturation 
scale), we assume a common tq for ^^snn = 130 and 200 
GeV. We note that final particle multiplicity is not sen- 
sitive to the initial time tq within a range 0.5 < tq < 1.0 
fm/c compared to the energy dependence of the number 
of initial gluons. 

Let us recall the relations among thermodynamical 
variables, i.e., energy density e, temperature T, and num- 



ber density n for the massless free parton system: 

^'-d.+d^mT^ 



2N,NsNf 



36, da 



2iN^ - 1) 



(10) 

(11) 
= 16, 



where dq 

C(3) = 1.20206, and S(= 486 MeV/fm^) is a bag con- 
stant which is fixed by matching pressure between the 
QGP phase and the hadron phase. Thus, we obtain the 
initial condition in the three-dimensional (3D) space from 
the number density (IC-71): 



T{to,x) 



TT^n{To,x) 



:S\ ^ 1/3 



43C(3) 

or, from the energy density (IC-e) 



(12) 



, 3Q[e{To,x)-B] \'/^ 



When the temperature at (tq, x) is below the crit- 
ical temperature Tc = 170 MeV, the all thermody- 
namic variables in the site are set to zero. Figure |21 
(a) shows the initial parton and energy distributions at 
the initial time tq = 0.6 fm/c as a function of space- 
time rapidity r]s. Here dN/dr/s ~ tq J d^x±n{TQ, x) and 
dEx/dris ~ Tq J d'^x_i_e{TQ, x). It should be noted that 
the assumption of the thermalization for gluons produced 
from the CGC is to reduce the transverse energy per par- 
ticle from Ex/Ng = 1.6 GeV to Er/Ng - 1 GeV within 
our present parameters. This change should be obtained 
by non-equilibrium descriptions of pre-thcrmalization 
stages. We also plot an initial condition which was em- 
ployed previously in Ref. js^]. Here we integrate the en- 
ergy density over all fluid elements with T > 100 MeV. 
The previous initial condition also reproduces the rapid- 
ity distribution of final charged hadrons. For longitudinal 
profiles, it consisted of two regions, i.e., boost invariant 
region near midrapidity I'qi.l < 2 and Gaussian shape in 
forward/backward rapidity regions. Here ry^ is a local 
space-time rapidity [S^ which is shifted from the global 
midrapidity according to the difference of the thickness 
of colliding nuclei at x±. The CGC initial condition has 
no boost invariant region due to the x dependence of un- 
integrated gluon distribution functions. It is interesting 
to see the transverse profile from the CGC and to com- 
pare with other parametrizations since transverse presure 
gradient causes radial flow in hydrodynamics. Three ini- 
tial conditions in the transverse direction are compared in 
Fig. 121(b) . In Ref. [s^l , transverse profile for the initial en- 
ergy density scales with the number density of the binary 
collisions (b. c. scaling). On the other hand, the flat pro- 
file with smearing near the ed ge o f a nucleus is employed 
in Ref. [sj. Note that, in Ref. |8ll|. the initial condition is 
designed for 6 = fm collisions and that cylindrical sym- 
metry is explicitly assumed. For initial energy density 



6 



2000 r -r^^ 2000 

1800- \^) dN/dris i1800 

1600^ dET^/dHs (flat+Gaussian) ^1600 

1400 
1200 



1200 S 




1400 ^ 



50 
45 
40 
_ 35 
E 30 
25 



> 
S20 



0) 



15 
10 
5 




-8 



L (b) 


— CGC 

b.c. scaling 


L b=2.4fm 


flat 





-2 2 
x(fm) 



FIG. 2; (a) Space-time rapidity dependence of the initial conditions dN/dr]s and dET/dr}s in Au + Au collisions at y^sjvjv = 200 
GeV. We also compare the initial energy density distribution with the one which we employed previously in Ref. |80|. (b) 
Comparison of the transverse profile from IC-n with those in Refs. Bind 



taken from Ref. |8l|, we take account of correction fac- 
tors 1.1 and (1.0/0.6)''/3 

coming from the differences of 
collision energy and initial time respectively. The value of 
maximum energy density varies from 13 to 43 GeV/fm'^. 
However, these initial conditions work at least for pion px 
spectrum by adjusting thermal freezeout temperature. 



III. HYDRODYNAMICS AND PARTON 
ENERGY LOSS 

In this section, we summarize essential features of the 
hydro-hjet model [H El H . 

Once an initial condition in a hydrodynamic simula- 
tion is assigned at an initial time tq, one describes the 
space-time evolution of thermodynamic variables by nu- 
merically solving hydrodynamic equations d^jT^'^ ~ 
with the help of the EoS, P = P{e). We are assum- 
ing baryon free fluids here, so pressure P is a function of 
energy density e only. In the ideal hydrodynamics, the 
energy momentum tensor is T^" = (e + PJu^u"^ — Pg^'^ , 
where is the four fluid velocity. Assuming mass- 
less free parton system as discussed in the previous sec- 
tion, the EoS for the QGP phase is P = (e - 4P)/3. 
For the hadron phase, we employ the partial chemical 
equilibrium model [s^l which describes the early chemi- 
cal freezeout picture by introducing chemical potential 
for each hadron. Both the phase transition tempera- 
ture and the chemical freezeout temperature are taken 
as Tc = T''^ = 170 MeV. As a usual prescription, the 
bag constant B = 486 MeV/fm'^ is chosen by matching 
condition for pressure Pqgp(Pc) = ^'hadron(Ic)- It should 
be emphasized that, at collider energies, relativistic co- 
ordinate in the time and longitudinal directions (t, ?/s) 
is a substantial choice to describe the evolution in the 
whole space-time as discussed in Ref. ^^Is^. Within hy- 



drodynamic approaches, the studies of centrality and ra- 
pidity dependences of observables simultaneously at the 
collider energies are only accomplished by fully 3D hy- 
drodynamic simulations with the relativistic coordinate. 
The statistical model and blast wave model analyses tell 
us that the kinetic (thermal) freezeout temperature T^^ 
is smaller than the chemical freezeout temperature T'^^ 
[S^ . Whereas, both thermal and chemical equilibrium 
is assumed in the conventional EoS. This means that 
both frcezeouts happen simultaneously and that, conse- 
quently, one cannot reproduce the px spectra and the 
particle ratio at the same time. By taking account of the 
chemical potential accompanied with the number conser- 
vation of each stable hadron between chemical freezeout 
and kinetic freezeout, we obtain the EoS denoting the 
picture of early chemical freezeout. Note that this EoS 
results in the almost T'^ independent pT slope for pions 




For the hard part of the model, we generate momen- 
tum spectra of hard partons by PYTHIA 6.2 [s^l in which 
2-^2 pQCD hard processes are included. Initial and 
final state radiations are taken into account for the en- 
hancement of higher-order contributions associated with 
multiple small-angle parton emission. EKS98 nuclear 
shadowing |84j is used assuming the impact parameter 
dependence as in Ref. js^. We employ the model in 
Ref. [s^l to take into account the multiple initial state 
scatterings, in which initial fc^ is broadened proportional 
to the number of scatterings. Initial transverse positions 
of jets at an impact parameter b are determined ran- 
domly according to the probability specified by the num- 
ber of binary collision distribution. Jets are freely prop- 
agated up to the initial time tq of hydrodynamic sim- 
ulations neglecting the possible interactions in the pre- 
thermalization stages. Jets are assumed to travel with a 
straight line trajectory in a time step. 
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Jets lose their energies through gluon emission induced 
by the dense medium. A lot of work has already done by 
many authors [H IH [Hil 113, lli ■ Here we employ 
the Gyulassy-Levai-Vitev (GLV) formula based on 
an opacity expansion approach which is a relevant for- 
malism for heavy ion collisions. In the original formula, 
the amount of energy loss is expanded into opacity. Here 
we take the first order term and neglect the kinematics 
of emitted gluons 

Ai? = C drp (r, X (r)) (r - tq) In (^^) , (W) 

where C includes the strong running coupling and color 
Casimir factors and p is a parton density which is to 
be taken from full 3D hydrodynamic simulations in our 
approach. /.j(=0.5GeV) is a screening mass and L(=3fm) 
is a typical length of a medium which can be identified 
with a lifetime of partonic phase, x and po are a position 
and an initial four momentum of a jet. The existence of 
T — To comes from the so-called LPM effect |83- In the 
static medium case, po • w is replaced by an initial energy 
of a jet Eq. 

In the actual simulations, the coefficient C in Eq. 1)14(1 
is regarded as an adjustable parameter roughly corre- 
sponding to an averaged value over quarks and gluons 
up to some factor. Our strategy here is that a constant 
C is set by fitting the PHENIX data on Raa for neutral 
pions in the most central events at midrapidity [T^ and 
that we next look at the centrality and rapidity depen- 
dences. The main reasons are the following. First, we 
neglect fiuctuations of the radiated gluon number [83 in 
our calculations. An absolute value of the energy loss is 
about factor two smaller when one takes into account the 
fluctuations of the radiated gluon number at RHIC [s^ ■ 
However, the pr dependence of Raa is found to be the 
same at RHIC energies [s^- Second, in our present hy- 
drodynamic calculations, it is assumed that massless free 
gas for quarks {N f = 3) and gluons. However, it is com- 
monly believed that gluons arc dominant components in 
the early stages of collisions at collider energies because 
of the gluon dominance of the parton structure function 
at small- a;. Thus inclusion of specie dependent factors in 
C may lead to discrepancy between our hydro-|-jet results 
and data. 



IV. RESULTS 

In this section, wc first study the centrality dependence 
of longitudinal and transverse distribution. Next, the re- 
sults for the centrality dependence of the nuclear mod- 
ification factors and back-to-back correlations are pre- 
sented. Finally, we study the nuclear modification factor 
at forward rapidity. 

Let us summarize our model parameters before show- 
ing the results. In the following, we use the parame- 
ters K = 0.7 and A = 0.2 in the gluon distribution of 
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FIG. 3: Rapidity and pseudorapidity distributions of all 
hadrons in Au -|- Au collisions at ^/snn = 200 GeV at 
h = 2.4 fm are compared to the corresponding transverse en- 
ergy distributions obtained from hydrodynamic simulations 
with CGC initial condition. 



a nucleon. k depends on how to match initial distribu- 
tions. For IC-e (IC-n), = 2.25(3.6). In the hydrody- 
namic calculation, the kinetic freezeout temperature of 
T'h = 100 McV and the initial time of tq = 0.6 fm/c are 
taken. We will not analyze the detailed T**^ dependence, 
but will make a short comment on that. A parameter 
C — 0.35 in the parton energy loss formula is fixed by 
fitting the nuclear modification factor for neutral pions 
from PHENIX data at 10% centrality 



A. Rapidity, centrality, and energy dependences of 
hadrons 

(Pseudo)rapidity distributions of all hadrons and 
transverse energy in central Au+Au collisions from the 
hydrodynamic simulation are shown in Fig. |31 Here 
IG-n is used. As one can compare Fig. 13 with Fig. [21 
(a), the effect of the hydrodynamic afterburner is to 
reduce the transverse energy per particle due to pdV 
work. We find that {dEr / dy) / {dN/ dy)\y^Q ^ 0.54 GcV 
and {dET/dri)/{dN/dvi)\^=o = 0.56 GeV. Recah that our 
CGC initial condition from Eq. (QJ yields Er/Ng = 1.6 
GeV and that it becomes Ex/Ng = 1.0 GeV after assum- 
ing a thermal state from Fig. |21 (a) . The consequence of 
the hydrodynamic evolution on the shape of rapidity dis- 
tribution is to make it wider very slightly. Final shape of 
the rapidity distribution looks like Gaussian which is con- 
sistent with recent data from BRAHMS (q^I- Our result 
supports that the KLN calculation jSJ which is based 
on the assumption of parton-hadron duality is a good 
approximation for the discussion of the shape of pseu- 
dorapidity distribution. The number of initial gluons at 
midrapidity is ~ 1480 (see Fig. 121(b)), which can be com- 
parable with the number of final hadrons at midrapidity 
(^ 1160). One see in Fig.Othat pseudorapidity distribu- 
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tions for both number and transverse energy are almost 
flat at mid-rapidity in the range of \r]\ < 2, although 
rapidity distribution has no flat region. 

The main difference between the previous and 
the present initial conditions is whether boost invariant 
region exists or not (see Fig. El (a)). One may ask why 
the previous initial conditions with boost invariant region 
near midrapidity can also result in the Gaussian shape 
in rapidity distributions, ft would be instructive here to 
discuss about the boost invariance for thermodynamic 
quantities and particle distributions within a hydrody- 
namic framework. Non-boost invariant dN/dy recently 
observed by BRAHMS ^ reminds us the so-called Lan- 
dau picture Boost invariance for dN/dy is equiv- 
alent to the boost invariance of thermodynamic quanti- 
ties in the original Bjorken assumption |6l| in which the 
length in the rjs direction is assumed to be infinite. How- 
ever, this is not always true if boost-invariant region for 
thermodynamic variables is finite. In such a case, rapid- 
ity distribution can have no boost-invariant region. Let 
us consider a simple example. Rapidity distribution can 
be related to the phase-space distribution function via 
Coopcr-Fryc formula |93 |: 



dN 
dy 



oc 



d^PT J f{x,p)p^'dafj 



(15) 



Consider, e.g., a finite Bjorken cylinder | rjs |< 770 with 
a radius R for massless pions. Assuming no transverse 
flow, we have / dcr'^ = / drisTTR'^T{coshr]s,0±,smhr]s) in 
the case of flat temperature profile. Thus the rapidity 
distribution becomes 

dN f°° 

— oc / drjs dpTPr cosh{T]s -y)f{p) 



-no 

(X / dr], 

-no 



cosh^(77s - y)' 



(16) 



where /(p) = l/(exp[pT cosh(77s — ?/)/T] — 1). This shows 
rapidity distribution for massless pions in a fluid ele- 
ment at rjs becomes cx l/cosh^(?; — 77s). After the above 
distribution is folded in the finite boost-invariant region 
I |< 7^0, one obtain 

^ cx T^[tanh{y + 770) - tanh(77 - 7;o)]. (17) 
dy 

The resultant rapidity distribution still looks like a Gaus- 
sian shape when rjQ < 2. The fiat region starts to appear 
when 770 >, 4. For more detailed discussions including fi- 
nite mass effects, see Ref. [q^- We actually reproduced 
the PHOBOS {^/^Wn =130 GeV) and the BRAHMS 
{y/sNN =200 GeV) data by u sing such a Bjorken like 
initial condition in Ref [13 El Eol (for an initial condi- 
tion, see also Fig. 121(a)). 

In Fig. ^ pseudorapidity distributions of charged 
hadrons in Au + Au collisions at both y^s^vtv = 130 
and 200 GeV are compared with the PHOBOS data 
Figures 4 (a) and (b) (Figures 4 (c) and (d)) show the re- 
sults of the initial condition which matches the number 



(energy) density distributions for partons between the 
CGC and hydrodynamics. Impact parameters in each 
panel are, from top to bottom, 2.4, 4.5, 6.3, 7.9, and 
9.1 fm (2.5, 4.4, 6.4, 7.9, and 9.1 fm) for = 200 

(130) GeV. These impact parameters are evaluated from 
the average number of par ticipants at each centrality es- 
timated by PHOBOS [9J|. We are not interested in fit 
for all rapidity region, because we do not include baryon 
chemical potentials and both the CGC and hydrodynam- 
ics break down when density becomes small. The results 
within the range I77I <,3-4 are satisfactory for reproduction 
of data. Therefore, KLN kr factorization approach in the 
CGC provides very good initial conditions for the hydro- 
dynamic simulations which reproduce rapidity, centrality 
and energy dependences. It should be also emphasized 
that it is not easy to parametrize such initial conditions 
which reproduce the data with the same quality as the 
CGC initial conditions presented here. 

In the calculations, we use A = 0.2 to get the best 
fit. However, we find that the results with the range of 
A = 0.2-0.3 are still within experimental error bars. Note 
that a good description of HERA data has been obtained 
by parametrization of the saturation scale ~ with 
A ~ 0.25-0.3 ISllsi. 



Kinetic freeze-out temperature T*^ which is conven- 
tionally fixed by the slope of the transverse momentum 
spectra is a free parameter in hydrodynamics. However, 
Pt slope for pions becomes insensitive to T'^ when one 
considers early chemical freezcout [s^jH^- But, we need 
to check the T*"^ dependence on the multiplicity. When 
we change the kinetic freeze-out temperature, we can also 
reproduce the data by = 3.5, 3.6, 3.7, and 3.88 for T*-^ 
= 80, 100, 120, and 140 MeV, respectively, for IC-n. We 
also investigate the tq dependence on the particle multi- 
plicity. We find that the particle yield with tq = 1.0 fm/c 
is only 4% larger than the case tq = 0.6 fm/c regardless 
of the centralitics as well as energies. This difference can 
be understood by the pdV work effect. 

Mean transverse momenta (pt) for pions, kaons and 
protons as a function of A'part arc compared with the 
PHENIX data 0| in Fig. Since the mean trans- 
verse momenta reflects the low pt physics especially ra- 
dial flow, contribution from hydrodynamic components 
at midrapidity is considered in this calculations. Al- 
though our results are slightly smaller than the data in 
central and semicentral regions, the overall trend is con- 
sistent with the data. For pions, the semihard spectrum 
starts to be comparable with the soft spectrum around 
Pt = 1.5-2 GeV/c in this approach. As we will see in the 
next subsection, we reproduce the px spectra for pions by 
including contribution from jets. So the deviation for pi- 
ons can be filled by the semihard spectrum. While semi- 
hard components for kaons and protons are very small 
in low and intermediate pt regions. A little more radial 
flow is needed to gain the mean transverse momentum in 
central collisions. More detail discussion will be given in 
the next subsection. 
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FIG. 4: Pseudorapidity distributions of charged hadrons in Au + Au collisions at y^sjvjv ~ 130 and 200 GeV are compared to 
the PHOBOS data [oj. For the choice of initial conditions and impact parameters, see text. 



B. Jet quenching 

In this subsection, wc compute the ccntrality and ra- 
pidity dependences of high pT spectra and nuclear mod- 
ification factors in the forward rapidity region. All re- 
sults are obtained by using the initial condition IC-n. In 
our model, high pT jets suffer interaction with the local 
parton density whose evolution is governed by hydrody- 
namics with the CGC initial conditions. We only take 
into account parton energy loss in the deconfined matter 

r>r,. 

Figure El shows that the pr spectra for identified 
hadrons from the CGC-f-hydrojet model are compared 
with the PHENIX data jOg for several centrality bins. 
The final pr spectrum in our approach is the sum of the 
contribution from fluid elements and from fragmentation 
of jets. First of all, we see that spectra obtained with the 
CGC initial condition yield the similar results with the 
initial condition employed previously |65l | . Agreement for 



the pion spectra with the PHENIX data is perfect. As 
mentioned before, hydrodynamic component is insensi- 
tive to the choice of T^^ for the pion spectrum in the 
chemically frozen hydrodynamics. Slopes of kaons and 
protons in the model prediction are steeper than that of 
the data. This indicates that more radial flow is needed 
for kaons and protons to reduce yields in the low px re- 
gion and to enhance in the intermediate px region si- 
multaneously, even after inclusion of the jet components 
into hydrodynamics. We have checked that freeze-out 
temperature of T^^ = 80 MeV still underestimates the 
slope. RHIC data suggest that initial transverse flow at 
the time of thermalization tq seems to be important as 
studied by Kolb and Rapp [93. This observation is nat- 
ural because one expects that many rescatterings lead 
to thermalization. Initial flow profile can be calculated 
from the classical Yang-Mills fields. Therefore, imple- 
mentation of Yang-Mills results into the initial condition 
in the hydrodynamic simulations could considerably im- 
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FIG. 5: Mean transverse momenta for pions, kaons and pro- 
tons as a function of A'^part. Here contribution only from hy- 
drodynamic components is taken into account. T'*' = 100 
MeV is used for all centralities. Data are taken from Ref. l9al 



prove (modify) the hydro results on the transverse mo- 
mentum distributions. Discrepancy in the high px for the 
proton spectrum might be improved by considering the 
recombination among hard partons and soft partons ■ 
In Fig. 13 we investigate the centrality dependence of 
the nuclear modification factor for \ri\ < 0.35 



Raa{pt > 4.5GeV/c) = 



4.5GcV/c 



dN 



AA 



5GcV/c 



dNPP 



(18) 



together with the PHENIX data |9||. Our results only 
account for the data up to mid-central events and fail to 
reproduce the data in peripheral collisions. Although this 
is reasonable because neither CGC nor hydrodynamic de- 
scription can be reliable in small particle densities, let us 
study more detailed mechanism. 

In the recent studies [Tsl IT^ l99l |. it is claimed that 
the centrality dependence of single and dihadron sup- 
pression spectra even in peripheral collisions is described 
by jet quenching in the partonic phase. They assume in 
their calculation that the parton density is scaled with 
the number of participants and that the system expands 
according to (l-l-l)-dimensional Bjorken expansion |6ll |. 
Line integral for the loss of parton energy is evaluated 
at all densities. This means that there exists only the 
partonic phase even in very small partonic density and 
that the effects of phase transition or the hadronic phase 
are neglected in those calculations. In order to study 
discrepancy between their results and our results in pe- 
ripheral collisions, we also estimate the centrality depen- 
dence of Raa assuming (l-l-l)-dimensional Bjorken ex- 
pansion with the same initial conditions as the 3D ones. 
In this calculation, the transverse shape of the initial con- 
dition is unchanged and the magnitude of the parton den- 
sity decreases as p{t,x±) = {tq/t)p{to,x±). Note that 



our transverse profile scales as p « Q'iib) /cts{Q'i{b)) w 
Ppa.it/ C(s{Qlib)) at mid-rapidity. pQCD calculations with 
(1+1)-D Bjorken expansion are performed until the max- 
imum partonic temperature reaches T = 100 MeV. When 
jet quenching happens in T > Tc, the result of this (1+1)- 
D Bjorken calculation yields almost the same results as 
our full 3-D calculations as plotted in the dotted line in 
Fig. \7\ However, when we do not impose a restriction 
on temperature in the parton energy loss as T > Tc, 
the result becomes closer to the data (see solid line in 
Fig. [Tj). The difference between these two calculations 
suggests that there exists a contribution to energy loss 
besides the purely partonic phase. This can be under- 
stood as follows. The density of the system essentially 
decreases like ~ 1/r even in our 3D hydrodynamics. In 
the case of most central collisions, initial density is so 
high and most energy loss is likely to occur in the early 
stages of the collisions. On the other hand, initial energy 
loss alone is not enough to lose energy of each parton 
due to small life time of the partonic phase in periph- 
eral collisions. Hence energy loss should last longer until 
the density of the system becomes small. From above 
considerations, hadronic interactions |lO0l | seem to be re- 
quired as well in order to understand the suppression in 
peripheral collisions, when the bulk matter scales as the 
dN/dy ^ Npaxt/ oisiQl) which is consistent with the cen- 
trality dependence of the experimental data contrast to 
the number of participants scaling. 

In Fig.lHl we plot the ratio of the strength of correlation 
in Au-f-Au and in pp collisions 



Iaa 



/f^d(A<^)Cf(A</.) 



(19) 



as a function of the number of participants A'^part together 
with the experimental data from STAR C2(A0) is 
the azimuthal pair distribution per trigger particle. 



C2(A(/') 



Atric 



gor 



dN 

I ^ dAcfidAi] ^ 



Here A(f) and A77 are, respectively, the relative az- 
imuthal angle and pseudorapidity between a trigger par- 
ticle and an associated particle. Charged hadrons in 4 < 

PT,triggcr < 6 GeV/c and in 2 < Pr.associate < PT,triggcr 

GeV/c are defined to be trigger particles and associated 
particles respectively. In the calculation of C2(A(/)), back 
ground is subtracted as in Ref. Near side correla- 

tion is defined as | Acf) |< 0.75 radian and away side as 
I A(j> |> 2.24 radian. Our results for the near side cor- 
relation show Iaa ~ 1 for all centralities, because inter- 
actions of jets after their hadronization are not included. 
However the data show Iaa < 1 for the near side in 
peripheral collisions. The results on the away side corre- 
lation for A^part/^ 100 deviate from the data, which is con- 
sistent with the behavior of Raa as discussed in Fig. [7\ 
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FIG. 6: Centrality dependences of the pr spectra for (a) pions, (b) kaons, and (c) protons obtained from the CGC+hydrojet 
model are compared with the PHENIX data Kinetic freeze-out temperature T**^ = 100 MeV is used in the calculations, 
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FIG. 7: The CGC+hydrojet result on nuclear modification 
factor Raa{pt > 4.5GeV/c) in Au + An collisions at y'sNN = 
200 is compared with the PHENIX data |9l. 



FIG. 8: The CGC+hydrojet result on Iaa in Au + Au col- 
lisions at y^sjvjv = 200 is compared with the STAR data [l^. 



C. Nuclear modification factors in the forward 
rapidity 

Recently, the pseudorapidity dependence of nuclear 
modification factors in both Au+Au and d+Au collisions 
observed by BRAHMS [Tsl Is^ attracts many theoret- 
ical interests in the context of the saturation physics, 
since one can go deeply into the small x region at for- 
ward/backward rapidity. Previously, we investigated the 
nuclear modification factor in Au+Au collisions 



RAAiv = 2.2) 
Raa{v = 0) 



(21) 



within the hydro+jet model |6J|. In Ref. |6J|, rapid- 
ity dependence of initial parton density was assumed to 
be flat in I rj'^ \< 2 and Rri was less than unity just 
because of the difference of pt slope of semihard spec- 
tra between midrapidity and forward/backward rapidity. 
However, the initial parton density from the CGC has 
no flat region. Therefore, we should revisit the rapidity 
dependence of the nuclear modification factor, because 
partons travel almost straight line and parton energy loss 
may directly probe the medium at each rapidity bin. 

The CGC+hydrojet result for charged hadrons is com- 
pared with the BRAHMS data [13 in Fig. EJ We ob- 
tain i?^ ~ 1 in the soft region pT < 2 GcV/c which is 
dominated by the hydrodynamic component. This is the 
same result as previous work even though initial rapid- 
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ity distribution is different. The reason is the follow- 
ing. Both initial conditions lead to reproduce the PHO- 
BOS/BRAHMS data on dN^h/di] in which multiplicity 
at 77 ~ 2 is almost the same as at 77 = 0. pr distributions 
at ?7 = and 2.2 are similar to each other below pT ^ i 
GeV/c in pp collisions [63|. We also find that the effect 
of radial flow at 77s ^ 2 is also the same as at midra- 
pidity. Hence Raa in the low px regions are insensi- 
tive to pseudorapidity within this range. One may worry 
about the effect of pseudorapidity dependent Jacobian 

J = p/E = \l ^ ~ JTi^/ (py cosh^ 77 -|- TTi^) on Pt distribu- 
tions. However, this is already canceled in the calculation 
of Raa before obtaining its ratio i?^ in Eq. H21|) . 

Situation is different in the high pT part, our result 
starts to deviate from unity at pT ^ 5 GeV/c since pt 
slope from pQCD hard collisions becomes much steeper 
for larger rapidity. On the other hand, experimental data 
deviates prior to our result. Note that the discrepa ncy 
between the present result and our previous result |6J| 
comes from the initial longitudinal shape of the energy 
density shown in Fig. [21 (a). In the previous parametriza- 
tion of the initial energy density, dynamical evolution at 
midrapidity was almost the same as at the space-time 
rapidity rjg ~ 2. Thus i?^ < 1 came from purely the 
difference of px slope between midrapidity and forward 
rapidity. While, in the initial condition from the CGC, 
initial energy density at midrapidity is slightly larger 
than at 77^ = 2. Therefore, the "slope" effect is com- 
pensated by the dynamical effect. This result indicates 
that Pt spectrum in forward rapidity region is already 
suppressed by the initial state effect. If this is the case, 
instead of using conventional pQCD approach, we may 
need to include x-evolution for the calculation of pt spec- 
trum up to Pt = 10 GeV/c that corresponds to x < 0.01 



at 7/ = 2 jsElzil- As demonstrated in Ref. study 
of the back-to-back correlations at forward rapidity will 
clarify the influence of parton saturation in the CGC. 



V. SUMMARY AND OUTLOOK 

We developed a dynamical model for the description 
of the time evolution of high dense matter created in 
heavy ion collisions at high energies on the basis of the 
CGC in colliding nuclei, hydrodynamic evolution for pro- 
duced bulk matter, and the parton energy loss in the 
medium. In this approach, collisions of two CGC's are 
calculated for the initial condition of hydrodynamic sim- 
ulations. Thermalized dense partons produced from col- 
liding two CGC's are evolved hydrodynamically and also 
used to evaluate energy loss of high-pT partons. 

We have shown that hydrodynamic simulations with 
the CGC initial conditions describe centrality, rapidity, 
and energy dependences of charged hadron multiplicity 
very well. The best agreement for the transverse mo- 
mentum distribution was obtained for pions with the hy- 
drodynamic component and the quenched hard jet com- 
ponent. However, our results on kaons and protons in- 
dicated that we need other mechanisms to understand 
kaon and proton transverse momentum spectra from low 
to high-pT- The centrality dependence of nuclear mod- 
ification factors and the back-to-back correlation were 
found to be consistent with the data up to the central- 
ity of A^part > 100-150. Finally we compared the nuclear 
modification factor at 77 = 2.2 to BRAHMS data by us- 
ing this new initial conditions. We found that, when the 
QGP parton density at the forward rapidity is smaller 
than at mid-rapidity, suppression of hard jet from pQCD 
calculation which is supplemented by DGLAP evolution 
is not enough to account for the data. It suggests that 
the small- a; evolution will be necessary for the calculation 
of particle spectra at forward rapidities at RHIC. 

Our main objective has been developing a consistent 
dynamical framework for the time evolution of both soft 
matter and hard jets. In order to meet this goal, many 
remains are to be done in the near future. Each compo- 
nent should be improved to obtain a better description 
of the whole of the stages of heavy ion collisions. 

For initial conditions, more sophisticated unintegrated 
giuon distribution should be used as mentioned in Sec.lTTl 
Numerical solutions of the classical Yang-Mills equa- 
tions in (2-1-1) dimension can be used for transverse 
profile in the boost invariant hydrodynamics. Quan- 
tum evolution is important for the LHC energy and for- 
ward/backward rapidity regions at RHIC. (3-|-l)-D Yang- 
Mills solutions in which renormalization group evolu- 
tion of the color charges is implemented are strongly de- 
manded. 

The ideal massless QGP EoS is employed in the decon- 
fined phase in our present hydrodynamic model. The EoS 
in hydrodynamics should be consistent with lattice QCD 
results. For example, one can parametrize recent lattice 
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results by a resonance gas model for hadronic phase |10 1| . 
In connection with jet quenching phenomena, it is argued 
that the small non-pcrturbative scre ening mass which is 
observed in lattice simulations |l02j leads to reduction 
of the amount of jet q uenching in the vicinity of phase 
transition region llO,*^ . Therefore, one expe ct th at jet 
quenching would strongly depend on the EoS |l04j . 

One of the important goals in the jet tomography is to 
constrain the evolving gluon density. Therefore, remov- 
ing the free parameter C in Eq. (|14|l will be one of the 
most important future work within our approach in order 
to see that gluon density provided by the hydrodynamic 
simulation is truly consistent with the energy loss of par- 
tons. In order to achieve this goal, we need to extend our 
hydrodynamics and use more complicated formula for the 
parton energy loss at the same time. For the more real- 
istic description of chemical compositions of quarks and 
gluons in the deconfinement phase, we can solve hydrody- 



namic equations together with rate equations for number 
densities of quarks and gluons [t^- This also affects jet 
quenching through the different cross section among gg, 
gq, and qq processes. The kinematic effect of radiated 
gluons, fluctuations, or higher order terms in the opacity 
expansion in GLV formula should be taken into account 
in the energy loss of hard partons. 
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